


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1988 


A comparison of two frequency domain 
adaptive beamforming algorithms for sonar 
Signal processing 


Nikitakos, Nikitas V. 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/22959 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist ae Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ies) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


" [ml ‘eta Py i if ry Py ry ee te ie RD ie) os be Sl le di tah i i ee | ale ete dd 
< ae Ma oo late ae a sag 6 8 Ng ol i. s4a0lu yp hee ee ee oh ey oy ; 
> riot ae vere h gictan HERAT AES a ee Se oe eee Se eee ery ere a ta 
*e ue Pt ates , ae one te) Ae Ce i Oe OOOO te Se ena nae eee eee ee ede uh yk ton 
Pr pa L A bt i ec eo oe ne we ei Te eerie ee Or Cit Seo ee fol aT YY) 

o Ch meme gel Si ; oe weeny Naan 2 aad Uy VL eed oo er ire tr Ie ON toa) 

' “ . 4 A A aed OCD) ew ct ee wy, a Ci eh. ei ie en i ey Te Se Ln Pele ty feel 

AL ! Cera ee fr ora Aye het 72 en OR Ot Ue Se) eee et ot tr 

ac , PACT ies Sa a) Sie ono i; FY ve ee ee © Rd Oe Aree FRPEee. ATE A we eee eo | iT) Yee tye yy tr 
A A i - ee cl eer yer} Ft ete TY aXe 5 Dn oe 4b FAL te 8 Steet trie eR TE eee ts 

0 7) rf ' an) thd ie PCa re | | y @ atta en & Be recor a west) 

oat a i ae Yor) ae en ae tak) 

Si a ad Pe Ie, ee oe kh oe ee et 

rer we) PAINS ULEt — POW 8 & & Benepe Dede Arartom ts dare 0908 

COM eet Ck Oe ae eT eee Morr rg avr l] ae ets 

D ht i S| et no ee ey A RT Th eT ee 












































70) or 2 er ee oe 
a OT TE ta ‘ A i ee be i | 
Bette 1s Om a azte 
be 








VT 
eT he ee a Ce Ce eee a a ew) * Cry = 
f D ORS Por tie ey tlie, Ba telern ie” Sah RR adler teint pets ipl. | feats d 
iat & Cs Te a ete to ee ey eer vy 
a be ay hi CL mere ee mot 20 pt erty we | 
en | 

















| ane) 


Ce TMA. AISI a eh en 


Pc i 
a OO et 
Aer) °e iy 


hy 
Pi 4 or | 
« ® Jat. er a or Th ee be 
Sratse 8 Boga gat 2! ao Oe te oa Soto wrth Te 
ee re i cr ee ie ere Teo ot 
4 U te Be ve 0.0 7 0 A ta meme A Se i OT Td eh adel | 
i ee Ok en Oe On Pe eo ee reer ort ee er. 
i Th A eee ee a PTT tt) Le on Lk aT Pry re Pre i) eet er ee eC eo 
oT a! ie er he oe ere SC orn Om TT oe aoe tt) ret) ee Pe erent wort ir iy ri wey 
O ‘40 args SL rsh ie Lt ee yhoo ee ee Ce eS tery he) vee ey Yer ty ed 
© Ba MP Bs pa pepe te ee , a ee ek a TK errr th bite! 

L 1 p Ger 4 Oe 308 (eee we Le ee dehy Pe 
mee @. Ca a Lv ee ee eee ee | 
a >. amRes Se bea OY] rd 
riye CAR ee eon ek ee Lek in) 

' Van a} LT een Len Le hk ee a oe ee) 
rd 1 ty W . A ORO) a a 2 re oe eee ee Or ee er ier) ear Sate ee as Pat ot ens 
a a Cot ce | Parry | b tee ’ + oe Pa CL Pe Teter ry ern tae ar eo PU ie ee ee he er ee TL 
. ra a Oe nn ae a err or a we eer ee i Se © Or Yee Pe Ye Co) ee Soren 
yl eee ty 4 Pe Bas ate Bgeh 0 var de te egg 8s BAS i VL od oh oe CR i he LL 
‘ ee oa) Un eee ee Os Ce Oe Oe ne Writ t a) en POPE ey er ores) oh Pe rs 
wa 


Sa wag Oe a, ea | tor eyed, Ant oe, 
U oa Pia acu) 4 Pre Br EI eal ee ee ee oe re eon Pr arc er ere a ne OO 
5 Neererrit) tare tert 


aT] x hee BS UY a ey Th A] Tare ee ee 
Ce e e beatae Se Pt ee Oe 0 Rete eA Oe Lie hh i beh lk ee he eh he 
Dee ee ee ek memere 


‘ a bh % GA ee er) et a ee en er Tire ote Boy) rin mtd 
| eed J 4 ‘ hd Meco MP ee ee TO Ce Sh ee eS ee eo err Sh ee hoe ea Lol 
ae ee et aT | Fama crs : Lea Th ° Pi cy ee oe, | 1 ay TER INBAP TTR GR Doge ry Cer Mi Oe Dat Beet ete ®e whee Mare 4 
ry i Lh ta ee en SPS Seer eS ee att hy Pa a adhe bs Renate haloes lion dar lees dart 
D SORA roe Dies ny. Snore mere pe hae cs Bree F oH OrMOR HON, 1 MoWasrags® tea i heed 
ad a) ee ae Ye ee?) ee ee ee en 
a ny voy wi A peek hk Bad Se be tol 
Bre gy oe bre LS ere Tt eee Ot ee ee, et 
ites ye Fe a Shs: <del. dan etapa ineeie tin dhn th eciabeatien ay dali 
Oe Se ay .? Co Ce i ee ee 
ca TT ee ee ee Se Let Rie ert) oe Pn? ee PT kh Le ore 
EU B WAL As teres 6 g@ DB? Bp RMMarqietag un tee ee CC ae we Sh 
1 "SO Reapace argh og ti bere ee eh eh ht ee re ee he oe Str 
w Vere espe eto se og CO ee ee ee) ee ee ee oe 
. a oF Oi a Fe erte Fae e: et Le ee ee Nn te E bndn tl Bibel ies 
or i eC eevee Ped PR te ee he ee 
or © % 9 ae t eh eh len De he a 
Sheet hh oh, leat 





ttt 
Lae 



































er 
e on ee i | 






































































‘ 

eee 

re ee tT ee ey 
be 




















fe sana core atn' hr"? 
ert ae fa Oe ee ee Oe Pr 
Oren Petit ete ok dh ee ee et Pen ok 
o° 0.8 8 Bet ee eee Ot ae ee eee 
TR rth ee eee rt ee erty ie bet . CL Lie eC ie te te 
@ite ee i eR te tt Et © TT ME Mile tedd en tel ie had eh) ST a tT 
ae e ten ae ti te eee hk tL LT et erty das 
rls st eaten ie hesitate be Lt Let | 
ey eg | ae Be he ee 
A ee 7 a ee he tt | 


ee ed he 
OE tte be te al 
tomer Bre Goaye 
Ch Se er ere 
te de ET hee el 
PL Tr ia me oe et ert S| 
n | eh hh ee Treen t 
oe an oT 2 %. ee fh ee ST eek er. ee en Ly f Te 2 toe OT ee SO mero 
Re . LT ee | ar Co. Pe Pe Te ey ee iy ee tee tt Le eee TC 4 yer 
F; < ' Carry A ar ar eee ae - a og ‘ er ’ PPC ee eee Mey RC Ps? ee eee ry HPS er wD Sakhe 20 Dorsewee @ soerioesors 5 ads 
are nr Te i ed ‘ Cera idee ee a bt eT ee 
MU mam ad one eo Fy res eT, a ra a iat tee) ee ee 
eee ery * OAT hh ee he te i ery Th 
Cr ee ee coenf 07 We 3 New Be tphRotttens grettie « ap 
ay © 0 ° Fn Th ee td Se, Td oe er Pee TT ae 
Ah Ca ee ee eT a ee oh kn eS 
eh ae F iD «4% 
i he he ee hoy 
Prarie. eked teem aoe 
Pere ee Piers SC tag ee i) a. 
vere. i nny | ee i ee eee oe 
a ae pow teremye fae coe gee eerie 
ee ee ee tora 
Peet o te cot wen- Pit: fis 5 
a .. =e en 2 aa 
ar ee nh - SE ae 
Pd ‘“y Sa er? Pe Teo 
er ee i? eh he ne eet 
oo ns a ey ose 
aa. sy Sei Tid (ey a Prey 
ay —_ Be hie et oe eee 
Th SF gto pte ee attaks a 
Fee a i ry 
* 

























re ta ie 
@Vvete Bt BOew Ac® 
fees 8 %& 
oe O 








eon. Th 
24a, Pras tea, 






U O a 
ee oe ht) r 
’ ee 0 
o- ' oy 
4 
a fo ce 
A Le | 


















Wrews, Cs 2 | Be ar 
et en oe | 





a ee | 
b a Pe Te 
















A a a TT 

































th 

a Cts a 

5 Ce a Per 

Or ee re ets 

ce i ey ay ee) 
y 
















” 
Ls 


te Ce nods 
rs Co 7 ree 
eow- atin P 4 eee tam 
Pa a Al od a a 
ty e td 
ms es o 
a Ce od eel 
: a = 
ora - 
a *. [i Pd 
. oF et M eens oe” dmyerer'ers 
5 D rr ee ee 
ae el - Denil = . Ce aw Te 
ed Cal al LS 7 an rd 
OO tell all add St Td Cd 
Pd a i ed ee 
oa al cad tt Maitre sO at ee cl al 
a) Per ed ~ Clade Caled 
Pee AM cack eee it a PO ee i Le ded) Pe ie ool Le 
ML aD ee ll od hl Ld Pd - ar a ay Ls ol ad adhe od 
Cd dee le Oe ott ae Part. Y ne? eh Fh ide iad he ae 
Pt iid ety © stew Orr lees gue Fe od Bad ot ed ad “ew bllnen Lens 
een ¥e. 6 » oe eee y arm astgen’ ee ee td Cateteviet pre 
o ore CP ad bel L Ce ie dal od ain tall ctl oc 
as ae} er od ° tes 
a Mew Coy 4tre f ¥ 
oy, of i oe ate Ce? «shy Peete Rd ahd te taal 
A Ale eae f al a 
Ce _- Le a OY | Pe le Le all ee ele EE te dal a bd 
G fe rare u i to ures cet Aree eat ee etre Cl I eal oti te adh 
wo. Oe, ee a ee Le | A Perey De ee ed Perr) oda dagll Oe eal | 
tig claret 8 ae Lt Loco Lt © ge thc ets eters 
ow, ee YC ee oe) ty = na Agee OC el i | 
Derren Me ere tee Pad © ORL Ot Var et 18 BSNS OEE 
Dd ek ied ee eS D FP we tees ot hw 2 O40 af Hee del adnate 
od aie we tw gate oe ed I Meet theted cel oath ohana thi tel ada 
Oh a oD rand ae oe et od er Ce ad at A al aaa 
one . Pe oe | . ee od oe od Re ed id Pe ad Ol al al Rated telat lt 
Ce ee cee hee ae ee Oe a ee Pe et ad 
CU de ar i es fod Re ede Med LT ee er ee oe ee oll aod 
Leer | Ce i a oe otha) CelP rerio’ are weeeneetg tg Boe why yd be a Y ae atta hel del ake LP bel ah eel 4 heed 
1e@ verre 6 4; 6 wWegere FA he of See rnd Cae ah al 
ew 38 Pr ed oe er) Y , Ce i eo es eo ed Po ti let okt Fe el tah tel ah ol al 
OP oe te Cl tt Ol el A A ode ed Ce i eM eet dod ol oll a! Ld ee oie 
Se, 68 set etree os ee) Ty ta 1) od hl 
Lot OL Bl tel eh Cf ee ee ed Ce ee Ned lt ea 
Ld a a : Ce al ee ae el ee Batt a 
were oe = ore orate 
ee age eee ee ee dali 
Ch ee ey CP ee i Le ad wriete OP ee 
ie ae a ed te ee eek Oe el al ee a del ee at ak at ee Lae cal ted bells Ded 
+ Cae Oe ae ee al es nd el Oe ae a ee Mo ee hal rere 
Pai heh) er i ee Le SE od OE Me Pa eT a alls 
re! Ferd OL ete pte ete ite arate ate ad 2 eT | eT Petite te 
ate ee CO he ee | STOTT Gg Glyde Bret ET ets WHET FRO H 1 FOLD Mate Pye! Coos *stowe¥en ¢ 
Pe Pr ee ee ee ee a ee 
Aa Pee ere eo ee nd ty 
Ot el dh hed Mite ca Kad atl octet tall ttc 
Py LD ee oe ed Pia wiewet Sal eelyp ganar arbi duchte, Aaa 
a ie ed nay one gtacaet 7 wees FS Hee.c are gs Cr Fd detetacth: wit naih el ak steal tnatel aati at ala 
'e *e #% Beer a] vee pre ng Ms oo 8 oT WP ee Oe ee od ot Pape arr pte acres cal 
14. ie Jes D oo r . Ss “ih weve @ OF LT ae ee et at a tO ee tome —_ WR treme 
. . 5 bint, SN OS SIGS ae oe On i. ee oe Me ee ee | i oe a PPP tt, ee a at) ee 
O Jonah ‘| a) ° Pe ee ee Te] Axe eb yne’ g Oe OFA Ve meh 18 ot # ot ae.e Oem O° ow ogy hiew OF a0 0 01 GeEe bt re caer nee ee 
¢e a : pe & ote a J CUS ot ee ee Lee Pee tee ee ee ee eo | Oe che Ad tel) Fe Deh ad hele altel ell al ee dialed ed 
u peed LU Lr ey . ae Ce ae Ce eee oe ‘nah ih ll ced what PT te tld of 
eh ee ee er eT Cel) ie LY i ee Te TT on nti [A CL Le tL TET olde tee dled MD a Le hae eh ahd Sela besetneadil 
eS ac, ae Pier es Teer ait ae ee eee Ee a al el ah eaiathmeebell 
J a ri @ae sue rcs a Fe "e Peed ol se remy Dd ele ek CL te al lll al ol ele, 
ovr ees $ ee Pe i Me ae 2 aenps ON diel a Pleat eee ae Std ial a Nk 9 een eed Li de tl dl titel 
‘ . td i ae) a ye ee ete ae 2 « Ci a oe ae eT a Nal ck Be ek Bl a ad dak ak ae atl tid 
MeO a tie ee nt do ee ek a we eet Mi ok ep att mith tlenl ane in ahh etl 
to ee ee ed ee er ee ot oe hd Pe ed il tell he el al oot Rt ed 
a or) pied gir hy Pal od Di hat ek ted ie oe die hd oie hh) el ol ad at a el ole ie Pd Ee Pele eh ee te ee Ld al 
So CU ee ae ee ee es ee OPT de tt On oO eee Le ee ll eel Md ol Dad tl 
ree We Jie, eR Niet aye Cee orn ceeuy aPC coy kd ee et Tse ee ol hn el oe A oO eed el Ra a es 
tC hes Fi re iy nT mr See eT Te of ot OO Do re A a al aes ahead Wh 
Lip Vt) Ce er eT Parr ee ty ae cece rr ee eh wert Ce Leiner 
ee Ce Wait pac) bevere COM HM F ears, | tae rie, eee ee) te he i bt ne Mn CP Oe eee od ehh) le) Pe oe Le 
ed Car aed oe en | PC a a ey a an et Tt eal at Me oT) i 2 te co ers Caden erkede p rarechi teadt 
Ea Cre Wer PO on 1s 2 PaO La Coc Sad, Racket ak “ai hed + SE Lik Sk ial ea ot id ah het hd hal 
Soe a Cer ea Pe Oe CS re Cr Be Te eh alt ee a Deel intitle Mika dS tdell dd, a 
e tdi Ce ee es aa ee Pe et a Lo inet OP a eh Stebewyras'e'y tee me Ce ee a a 
i) Y= e Pai typ ote Sere Po Pte eed ee ee of te et eed ee a 
Ce at Me ee) St as Lee ee ct Soe | eo NII ee ee Oe CP te be te et Ie ee hd Cee i ea Pe te lara 
(hg ek A Be eh TA Cs a) FF yte eet ae 40h, ver he poe Shah ge  Aacre Mestensbrg es othe fs! Hon Mb hdhel-ad cheated saiabantiog 
en es Par Se i ee ee rs 2 Le Ty sae iat 1 ite o ohare F Yee Prd Re we te tate Po ee | ed ee RT ol 
oe oa Pt oe id rd 2 ate aus @ it a od hae as ed yd Ce Le i ld Cee DL eae 
See tho sax Pe Cite ee eee eT) ey MT ho ee eT ia oe be ee a ol 
A ne, ire Sc i a nt en ore) Omen Cr tes Tete ot nm Me Okey 00 nen ene akin od 
Cn ' Bae feo a) shor 88 Le eve LY Pe ls eet ae fe te eee Oe ee eed a eee as 
ry ‘ en Me er 2 ee ee Tere roe Tia irr Ee eh Peete etl eee te eee eee ee ee 
ahs Oe en eter Po POV Fete eter“ OT A ae ee ee gee ait Pree 68 At eertyerger opr pt a ems 
CP eee Oe ee er Te TY ee oe ee i te te ee er er 
ee OF BIheyt tg te yey a See eid bt ae oe edd at tel) A ee Me he a! de Sh tat eee et Pte Le eas 6d Td 
r | ee i eZ » bare ae SY oY or em ad CA et em RY Ae eee ade th Oe, ei te ek 
CA i ie TO Le a a A ae he Bd i et ek kT A dk Sd a ee eT erat) 
Cr ee Le Te Ok ee ie tee Mt Fe ey ymrrye «Pattee gelewt sg ace @ gy prese ht hoped beg fg “er aa a paca amt pit 
U eset 8 UT A en MT eh sO he SA Oe el kD ee th A ahh etl Ri Lied al 
o a} ° Ft Mohd ee 4% am AT ea POUT SS or ok ee ee ee aul Dk De Ae ala cadet had tented 
CO Pr eS oe ae ade he a se Ne a ee ad | cece Tee eee ree EO TT os tees 
i ee Co Pre ee eee a) Thee ht. ee at ee 2 te tee rari et Oe a ae De hol a 
ers ee he Ce Se be Ce hy a far 2 OEM Ly LL) 20 te 21122 I LA RA oi le Dek Se ke cof lec atta ah cee bi dh eas he Fs, en 
a ee ot fe tatet Err rt fae ernie Meer 1) Dien Ea aks te cell ke be Aen AT lial af ee 
re a : Pere EMEP eT Tet at the hd eke an oe Lt ei A Yd Spe end pr 
’ rt ody i and Ln CL A oe Bi hth biel titel sie ad eee Ma eke 
ta s PPC irr os Ate ROTH LG EPP eee FTL Oey eos W ree gIMR pee adOr TD 
Le OL es It ae Pied terisy ach ok eat Le ld banat be ek I a etal ed os eet 
rey Ne i ie ir teee tet teria iar te gk edad Ue fod ed I ee ee od Ta eo 
CP A Joe Peer te ar ee Pr PESO Le ee ee Red ad RD oe oe a aha ae ali ok tn eared 


































































Ce | 
ed ors ore 
Ch ae er teuegnraee Pees ee Oe ed 





OO eo led ed ol ol ed A 









C.  ) on Coe Si Dal Gy hie 













































































vit ie aE) 















$ ia ? ose ete eg at epee ee ee er rrreeey ie Yer ne yee Cok ee ee ad ead read 
. i a er ie) “ ri ahh ay Pave FF wteer ath Aa wee See ays omer Depo sate Pad Gh Wk Be bd Se ol ae oa Oe eat te te 
° a) ’ A ae] i ree moa) = ee Oe ea te Lr) ee er ee Ree Te i PR ee An ee a a ee ail 
o ° -¥ y Pri ne a ers Ce MP eT Te toe ated RO 0 PR NL oll ot al AOR Sh dated alte dapat Sy PSO bcs dapat 
e a aed a 7 I fl gests 6G PW Me PT AT Lee la Oe On Wie or a) fe OL ok ey se ppt adit a! bit Deh te tas 
- ICL UNM R HN ML MUNN I tA RNAS. Ctr ian atl eta Bedt nd rien eaten Th eo 

. Of ry ’ c °, oF hs ad 4 Ce lhe aie a eI i 
: iu SF | P aan ft en cot J Tk, Td lee Lk i at I? Mh a el fi adc iba lf cede a Be lb a eld tg 








oUJe 


iA 


: CAtLae U 














NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


A COMPARISON OF TWO FREQUENCY DOMAIN 
ADAPTIVE BEAMFORMING 
ALGORITHMS FOR SONAR SIGNAL PROCESSING 





by 
Nikitas VY. Nikitakos 


a @ @ 


December 1988 


Thesis Advisor Lawrence J. Ziomek 


Approved for public release; distribution is unlimited. 


S15 Thy / 





Teil) fa V 


nclassified 
curity classification of this page 


REPORT DOCUMENTATION PAGE 


1 Report Security Classification Unclassified 1b Restrictive Markings 





1 Security Classification Authority 3 Distribution Availability of Report 

> Declassification Downgrading Schedule Approved for public release; distribution 1s unlimited. 
Performing Organization Report Number(s) 
a Name of Performing Organization 7a Name of Monitoring Organization 

-aval Postgraduate School (if applicable) 62Zm Naval Postgraduate School 





: Address (city, state, and ZIP code) 7b Address (city, state, and ZIP code) 

fonterey, CA 93943-5000 Monterey, CA 93943-5000 

1 Name of Funding Sponsoring Organization {8b Office Symbol 9 Procurement Instrument Identification Number 
(if applicable) 


> Address (cin, state, and ZIP code) 10 Source of Funding Numbers 
Program Element No Work Unit Accession No 
| Title (include securify classification) 


\ COMPARISON OF TWO FREQUENCY DOMAIN ADAPTIVE BEAMFORMING 
ieee VS FOR SONAR SIGNAL PROCESSING 
? Personal Author(s) Nikitas V. Nikitakos 


3a Type of Report 13b Time Covered 14 Date of Report (year, month, day} 15 Page Count 
faster’s Thesis From To December 1988 it 


} Supplementary Notation The views expressed in this thesis are those of the author and do not reflect the official policy or po- 
tion of the Department of Defense or the U.S. Government. 


’ Cosati Codes 18 Subject Terms (continue on reverse if necessary and identify bv block number) 
‘eld Frequency Domain Adaptive Beamforming Algonthms-Planar Arrays -Sonar signal proc- 





? Abstract (continue on reverse if necessary and identify by block number) 

Computer simulation studies of two frequency domain adaptive beamforming algonthms are presented. The two algo- 
thms are the frequency domain adaptive beamforming modified least-mean-square alzgonthm, and the proposed new fre- 
uency domain adaptive beamforming algonthm based on modified adaptive linear prediction-error filtering. The simulation 
udies were conducted to determine the multiple broadband target localization capability and the full angular coverage ca- 
ability of the two alzgonthms. The number of iterations that the adaptive algonthms took to reach a minimum estimation 
ror was determined. The algorithms were evaluated at several signal-to-noise ratios. Finally, using the results of the simu- 
ition studies a comparison between the two algorithms 1s performed. 


) Distribution Availability of Abstract 21 Abstract Security Classification 

J unclassified unlimited [} same as report CL} DTIC users Unclassified 

2a Name of Responsible Individual 22b Telephone (include Area code) 22c Office Symbol 

awrence J. Ziomek (408) 646-3206 


Gemeey 1473,84 MAR 83 APR edition may be used until exhausted security classification of this page 
All other editions are obsolete 





Unclassified 


pe - 


Approved for public release; distribution is unlimited. 


A Comparison of Two Frequency Domain Adaptive Beamforming 
Algorithms for Sonar Signal Processing 


by 


Nikitas V. Nikitakos 
Lieutenant, Hellenic Navy 
B.S., Hellenic Naval Academy, 1980 


M.S. in Applied Mathematics, NPS, June 1988 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCILNCE IN ELECTRIGAL ENGIN EEIING 
from the 


NAVAL POSTGRADUATE SCHOOL 
December 1988 


ABSTRACT 


Computer simulation studies of two frequency domain adaptive beamforming algo- 
rithms are presented. The two algornthms are the frequency domain adaptive 
beamforming modified least-mean-square algorithm, and the proposed new frequency 
domain adaptive beamforming algorithm based on modified adaptive linear prediction- 
error filtering. The simulation studies were conducted to determine the multiple 
broadband target localization capability and the full angular coverage capability of the 
two algorithms. The number of iterations that the adaptive algorithms took to reach a 
minimum estimation error was determined. The algornthms were evaluated at several 
signal-to-noise ratios. Finally, using the results of the simulation studies a comparison 


between the two algorithms 1s performed. 
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I. INTRODUCTION 


One of the major research items in underwater acoustic signal processing is to de- 
velop new sonar signal processing algorithms capable of quickly and accurately solving 
target localization problems. Present technology dictates that several lines of bearing to 
a target must be obtained before a sonar fire control solution can be computed. Ob- 
taining these lines of bearing is a time consuming and often dangerous task due to the 
increased probability of counterdetection and, as a result, evasive maneuvering and de- 
fensive action on the part of the target. A sonar system capable of providing timely, 
accurate target localization while minimizing own ship maneuvering would result in 
longer firing ranges and, therefore, a reduction of the threat to one’s own ship. 

Frequency domain beamforming is accomplished by applying appropriate phase 
shifts at the sensor outputs of an array to account for the relative propagation delays 
of a signal from a particular direction. The phase-shifted signals from all sensors are then 
added together coherently to realize the full array gain. Discrete Fourier Transform 
(DFT) beamforming is the usual method of determining the direction of arrival of a 
plane wave signal. A discrete number of direction cosine bins are formed and each bin 
corresponds to a discrete direction. If the number of direction cosine bins 1s large. very 
fine spatial resolution can be obtained. The phase shifts needed to cancel the relative 
propagation delays can be determined adaptively. 

The research performed in this thesis is a continuation of the work of Ziomek and 
Behrle [Ref. 1] and fully evaluates the capabilities of the complex, modified least-mean 
Square, frequency domain, adaptive beamforming algorithm that they developed via 
computer simulation studies. Also, a new frequency domain adaptive beamforming al- 
gorithm based on modified adaptive linear prediction-error filtering is proposed, which 
pronuses better resolution in low signal-to-noise ratio environments especially when 
multiple targets share common spectral lines [Ref. 2,3]. The computer simulation studies 
Were designed to test the two algorithms specifically for their multiple broadband target 
localization capability, their full angular coverage capability, and their angular resol- 
ution as a function of the input signal-to-noise ratio (SNR) at a single element in the 
array. Using lower SNR values than those used by Ziomek and Behrle {Ref. 1] in the 


evaluation of the modified LMS algorithm, and changing some assumptions. such as the 


number of array elements, sampling rate, and number of iterations, a comparison was 
performed between the two algorithms. 

Each target was modeled as a broadband sound source. As a result, the frequency 
spectrum of the output signal from each element of the planar array contains several 
frequency components. An estimate of the bearing and depression angles associated with 
each frequency component is provided as a result of processing the output frequency 
spectrum from each element in the arrav through the frequency domain adaptive 
beamforming algorithms. Therefore, if each target exhibits at least one unique frequency 
component (or spectral line), then all targets can be located. 

Full angular coverage is the ability to localize a target regardless of its relative po- 
sition with respect to the array. The broadside case is the easiest since it is at this posi- 
tion that the far-field beam pattern beamwidth is its narrowest. The endfire case is the 
most difficult since the far-field beam pattern beamwidth 1s the broadest at this point. 
The full angular coverage and multiple broadband target capabilities were tested simul- 
taneously. 

Baseline results Were the first assembled. They are defined as the bearing and de- 
pression angle estimation errors (measured in degrees) as a function of harmonic num- 
ber, and the number of iterations of the two algorithms, for the “no noise” case. 
Identical cases were then run using additive, zero mean, white, gaussian noise to corrupt 
the output signal from each element of the planar array. Average bearing and depression 
angle estimation errors Were plotted as a function of the input SNR at a single element 
of the planar arrav, and the harmonic number. 

Chapter II describes the theory used in the development of the two frequency do- 
main adaptive beamforming algorithms. First the already known modified complex 
least-mean-square (LMS) adaptive algorithm is summarized. Then the new algorithm 
which 1s based on modified complex adaptive linear prediction-error filtering 1s discussed. 
For both algorithms the construction of the direction cosine estimates, and the angle 
estimates will be presented in detail. 

Chapter III contains computer simulation results and an explanation of these re- 
sults. The results are presented as plots of the average estimation errors of the bearing 
and depression angles at three distinct signal-to-noise ratio levels (-3dB, -6dB, -9dB) for 
five different cases for each of the two algorithms. The different test cases include targets 
evaluated at broadside, endfire, random placement, and targets which share a spectral 


line. Conclusions concerning the effect of harmonic number, sampling rate, number of 


arrav elements, number of iterations used in the two algorithms, and SNR are made. The 
Appendix contains tabular numerical data for all results. 


Chapter IV will present final conclusions and recommendations for future research. 


II. ANALYSIS 


This chapter will present a brief analvsis of the two frequency domain adaptive 
beamforming algorithms for planar arrays which are evaluated in this thesis. The goal 
of both algorithms 1s to localize multiple broadband targets by processing the output 
electrical signals from a planar array of sensors to provide estimates of direction (both 
bearing and depression angles) and frequency content of the multiple acoustic fields in- 
cident upon the array. The development of the modified complex least-mean-square 
adaptive beamforming algorithm for planar arrays is a summary of the analysis section 
of a paper by Ziomek and Behrle [Ref. I]. However the modified complex adaptive lin- 
ear prediction-error filter beamforming algorithm for planar arrays is presented for the 
first time. As a theoretical background for both algorithms, the principles of adaptive 
beamforming for planar arrays in the frequency domain are presented in the next sec- 


tion. 


A. FREQUENCY DOMAIN ADAPTIVE BEAMFORMING FOR PLANAR 
ARRAYS 


Consider a planar array of M x N (odd), equally spaced, point source elements lying 
in the NY plane where M and N are the total odd number of elements in the X and Y 
directions, respectively. Let the random output electrical signal received at time instant 


fandelement ‘m,n, in the array be given by 


renin) = Vile) + nn). a — 1 0 
= ot Oe (2.1) 
ii=— Ne eOae 


Where y/dm,n) and nflm,n) are the deterministic signal and random receiver noise, 


respectively, 


ae (22) 


Me SO We: (2.3) 


ae (2.4) 


and 


Ker | (2.5) 


is the total number of time samples that must be taken per element in order to avoid 
aliasing when the signal y//,m,n) 1s composed of K harmonics [Ref. 4]. If we designate 
the length of the data record (1.e., the fundamental period ) recorded at each element in 
the array as 7, seconds, then the fundamental frequency (1.e., the FFT bin spacing ) is 
h= a Hz, the highest frequency component that can be calculated in y/lmnj is Kf 


I 
Hz, and the output received signals r/lm,n,; must be sampled at a rate of 


f, = L/T) samples’sec, (2.6) 


me@erest is given by Eq. (2.5) [Ref. 4]. Taking the DFT of Eq. (2.1) with respect to the 


time index / yields the complex frequency domain samples. 


Rigen) = Yigain) + N(@,Fitn), g=—L’,....0,...,L' 
m=—hKT....,0,...37 (257) 
n=—N",...,0,....N’ , 


Where the index q represents the harmonic number. If the acoustic field incident upon 
Y ¢ Py 


C 
rection (see Fig. 1), where 7, is a unit vector and g(t) is an arbitrary baseband function, 





the planar arrray 1s a single general plane wave g(i + ) propagating in the +”, di- 
then it can be shown that the output electrical signal at time instant / and element (m,n) 


in the array is given by [Ref. 4 p. 160]: 


uid, + vond,, 


vlan) = g(IT, + C 


Figure 1. 
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The corresponding frequency spectrum 1s given by [Ref. 2] 


J2ngfougndy _ p2rqforondy 
c 


(q.m,n) = Le, exp( + O(a ee (2.9) 


where 


Co = Y(qg,0,0)/L, g=— L’,...,0,...,L° (2.10) 


are the complex Fourier series coefficients that can be used to represent g(t) by a finite 


Fourier series with K harmonics during the time interval [7] < > h= 7 
0 


Ug = SIN Oy COS Wo ea) 


and 


Vo = Sin 8p sin Wo (le 


are the dimensionless direction cosines with respect to the X and Y axes, respectively; 
Q, and , are the depression and bearing angles, respectively, d, and d, are the interele- 
ment spacings in meters in the X and Y directions, respectively; and c is the speed of 
sound in meters per second. 

Equation (2.9) 1s the basis for the frequency domain adaptive beamforming algo- 
rithms to be presented in this thesis. The algorithms presented use a planar array to 
estimate both @, and y, for each harmonic q present in the multiple incindent plane-wave 
fields via frequency domain adaptive beamforming. Therefore, we shall use the following 


generalized version of Eq. (2.8) 


= Up,mnd + Yop, P 
ylLmn) = > eillT, a | (2518) 
kK 


which represents the output electrical signal at time instant / and element (m,n) in the 
planar array due to several incident plane-wave fields, where g,(7), u,, and v, are the ar- 
bitrary time function and direction cosines associated with the kth sound source (target), 


respectively. 


B. FREQUENCY DOMAIN MODIFIED COMPLEX LMS ADAPTIVE 
BEAMFORMING ALGORITHM 


The frequency domain adaptive beamforming algorithm to be discussed is based 
on processing the output complex frequency domain data R(q.m,n) from all M x N ele- 


ments in the planar array. We begin by defining the complex estimation error as 


e(q) = s(q) — $(q) (2.14) 
where 
M' N 
=a > D, IRgsmn) expl +i R(q,0.0)] (2.15) 


m=—Nf' n=—.N 


is the reference signal; 


At Nv 
n l 
=~ dD, Dd, amdlqnR(qunn) 
Sr 2S 4 
iE ) (2.16) 
c (q)R(q)d(q) 
Ley 


is the estimate of s(q); c(q.m) and d(q,n) are the unit magnitude complex weights in the 
X and Y directions, respectively; c(g) and dq) are the M x 1 and N x | complex weight 


vectors in the X and Y directions, respectively; given by 


c(q) =[c(g, — M’),...,c(g.0),....C(q, WM’) (2.17) 


and 


d(q) = [d(q, — N’),.-,d(q,0),..4(QN))" (2.18) 


respectively; and 
PQa we pai Jag, —O%,,0)...R(g, — MN ) 


R(G.0,0) —- —-R(q,0,4”) 


R(qg)=|  R(q,0.— N’) 
. R(qg,M',0) ... R(g,M',N’) 


Rigevie NO) 


is the M x N complex data matrix. Next, define the (M+ N) x 1 complex weight vector 


u(q) as follows: 





x(q) -| a | (2.20) 
Therefore. 
c(q)=A4 w(q) (2.21) 
where 
ee | a J (2.22) 
vif MaxN 


isa Mx(M+N) matrix, J is a M x M identity matrix, andQ is a Mx N null matrix 


as indicated; and 


aq) = B w(q) (2) 
where 
B=[ 9 | I J (2.24) 
NxM NxN : 


isa N x (M+N) matrix, Q is a N x M null matrix, and J is a N x N identity matrix as 


indicated. Substituting Eqs. (2.21) and (2.23) into Eq. (2.16) yields 


_ w'(QZgu(q) 


$(q) = “T1375 7 


where 
Z(q) = A’ R(@)B (2.26) 


is a(M+N) x (M+N) complex matrix. The complex weight vector that minimizes the 


mean-square error E{\e(q)|?} 1s given by [Ref. 1] 


Ww. 4(9) = wg) + 2uje(QLZ(q) + Z"(q@] w;, ¢ =0,1,2..... (22 
where 


e(q) = s(q)— $(q)_—- (2.28) 


is the estimation error after the ith iteration, 


wy (qg)Z(q)w{q) 


Gas (2.29) 


1s the estimate of s(q) after the ith iteration, and, the step-size parameter y,; 1s given by 


py ee te ee ee (2.30) 


Where yy is constant and o? and o? are the signal and noise power, respectively, at the 
center element in the array. After each iteration, each component of the complex weight 
vector w_,(g) 1s normalized by its respective magnitude in order to maintain unit mag- 
nitude, 

Once the complex weight vector w,_,(g) converges to a steady-state value w,,(q), the 
Steady-state complex weight vectors ¢,,(q) and d;.(g) can be obtained from Eqs. (2.21) 
and (2.23), respectively. The estimates of @,(g) and W,(q) at each harmonic q are given 
bx [Retas) 


éo(q) = sin f (ah igy + S(@ry?], q #0 (2.31) 


and 





nat f° (q) 
Wo(q) = tan AE5(a) g #0 (2.32) 
0 


where ui5(g) and 145(g) represent estimates of the direction cosines obtained by using a 
least-squares fit to the “unwrapped” steady-state phase weights 6§.(g,m) and $&,(q.n), 


respectively. Note that in the absence of noise 


27 Ofyito(q)mnd 
U o'0 x : ; 
O5s(qt) =e Le Opa @128) 
fe —ViewO on, 

and 

27afovo(q)nd 

U oo ; 

bss(qt) =) HL, 0, (2.34) 


The steady-state phase weights need to be “unwrapped”, that is, allowed to take on val- 
ues outside the closed interval [ — z, z], in order to ensure full angular coverage (1.e., 
0 < 6,(q) ay and 0< y,(g) < 2z) and correct depression and bearing angle estimates, 
6,(q) and w,(q), respectively. Finally, substituting Eqs. (2.15) and (2.16) into Eq. (2.14) 


yields the following expression for the steady-state estimation error: 


M N’ 
] 
ess) = Thy s yy |R(q,m,n)| exp[ +/2.R(9,0,0)] 


m=- MM’ n=—N’ 


— exp{ + j[05s(q.m) + b55(q,n) + LR(g,m,n)]} 


(2.35) 


where 


R(q,m,n) = |R(g,m,n)| exp[ +/2 R(q,m,n)] (2.36) 


css(qum) = as s(qun) expL +/055(qum)] = expl +/8s5(qvn)] (2.37) 
and | 

dss(qsn) = bss(q,n) exp +/Pss(qn)] = expl +ibss(4,7)3 (2.38) 
Where a,;.(q,m)= 1 and b,,(q,”) = 1 are real, unit magnitude, amplitude weights and 


6% (qm) and @%(qg,n) are real, “wrapped” phase weights. The procedure for the unwrap- 
ping of the steady-state phase weights in order to obtain least-squares estimates of the 


direction cosines at each harmonic is described by Behrle [Ref. 5}. 


C. FREQUENCY DOMAIN MODIFIED COMPLEX ADAPTIVE LINEAR 
PREDICTION-ERRROR FILTER BEAMFORMING ALGORITHM 


Consider a planar arrav of M x N (odd), equally spaced, point source elements lying 
in the AY plane where M and N are the total odd number of elements in the X and Y 
directions, respectively. Figure 2 shows a linear prediction-error filtering structure in the 
frequency domain (Le., after the operation of a time domain DFT). The frequency do- 
main adaptive algorithm to be presented is based on processing the output complex 
frequency domain data R(q,m,n) given by Eq. (2-7). Since the target localization prob- 
lem requires only phase informations 1n our algorithm, we processes phase only data. 

Let the one-step forward linear prediction errors for harmonic number q in the X 


and Y directions, be given by 


Fy (gst) = 5x(q,1,") — $,(q,in), nm=—N",...,0,...,N" (2:53)) 


and 


F, (qm,i) = 5(q,m,i) — §,(q,m,)),  m=— M’,....0,..,Af (2.40) 
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Figure 2. Frequency domain linear prediction-error filter structure 
respectively, with 
N= = (2.41) 
and 
M' = et | (2.42) 


where i is the number of the reference element in the direction of processing; s,(+) and 
5,(«) are the reference signal and the estimate of the reference signal in the X direction, 
respectively; and s,(+) and $,(+) are the reference signal and the estimate of the reference 


signal in the Y direction, respectively. Each one of the above signals is defined as follows: 


S,(q,i,n) = R(q,tn), (2.43) 


P 
Magi) => > AalqeR(qui + pr) =e Aalqurde(qisn (2.44) 
p=! 


S,(q.,i) = R(q,m,i), (2.45) 


and 


P 
| 
S(qn,i) = > ds (qum.p)R(qum,i + p) = = by (qm)x(qum,i) (2.46) 


where /,(g,n) and /i,(q,m) are the P x 1 complex weight vectors in the X and Y directions, 
respectively , and r(g,i,2) and r(qg,m,i) are the P x 1 complex frequency domain data vec- 


tors in the X and Y directions, respectively, given by 


Ag(g.n) = CAs(Qs1.2), Bel G.22) eer Pay) (2.47) 
Koa) = (Rogar Geet ree Rig,it P ny)’, (2.48) 
A,(qam) = (A, (qu). Ay(qutt.2)enee AGP), (2.49) 
and 
r(q,m.i) = [R(qym,i + 1),R(qum,i + 2),.......R(g, mi + P))" (2.50) 


Where P 1s the order of the linear prediction-error filter. 
Substituting Eqs. (2.43) and (2.44) into Eq. (2-39) and, substituting Eqs. (2.45) and 
(2.46) into Eq. (2.40) vields 


P 
Fy (qin) = R(q.i,n) — 7 > Alqpr)R(qi +p), n=—N’,...,0,...,N 5 ean 


p=1 


and 


P 
. l ° f ae 
Fy (qu!) = R(g,m,) — 7a > A(qm.p)R(qun +p) m=— M',...,0,..,.0. (2.52) 
p= 


Since we are working with phase only data, Eqs. (2.51) and (2.52) become: 


F,(q,in) = expl LA(G,i,")] 


EP 
| ae 
> > We(q.p.n)| expl +/LAg(9.p.0)] expl HL R(qi + pin), (229) 
p= 
ic — = ash Bao 


and 


F,(qm.i) =exp[ +/2 R(q.m,i)] 


iE 
7 > > Iy(qump)) expl +/LA,(q.m.p)] exp[ +/LR(q,m,i + p)]. (2.54) 
p=l 
m= —M..,.,0,...,-1'. 


From Eqs. (2.53) and (2.54) it can be seen that if each component of the complex 


weight vectors have unit magnitude, then the forward linear prediction error become 


zero when 
P 
exp[ +L R(q.i.n)] ae = » exp[ +jLh,(q.p.n)] expl +jL R(q,i =~ De). (2 55) 
PSL | 
i= —N",...,0,...," 
and 
P 
exp[ +/ZL R(g.m,i)] = 7 = exp[ +jLh,(q,m,p)] exp[ +/2R(q,m,i + p)), (2.56) 


AE AM. O00, -0 


Equations (2.55) and (2.56) will be valid only if the phase weights cophase the data at 
the output of each element in the prediction-error filters. 

The algorithm used for the optimum estimation of the complex weight vectors 1s a 
modified, complex least mean-square (LMS) algorithm. The estimates of the complex 
Weight vectors in the X and Y directions after the k-th iteration in X and Y directions 


is given by 


LA(9.) Ina oa [A,(9.n)]x =e 2u,LF, (qin)Iele (q.in)], 
== NY, 0)...0 (2.57) 
ke 0,12 330. 


and 


A, (q.m) car = Ly (qm), + 2uglF, (qn ale (q.,/)], 
m= —M'....,.0u,1" (2.58) 
A =0,1,2.3..... 


respectively. where the index k denotes the k-th step of iteration, F, (+), and. i -), are 
the forward linear prediction errors in the X and Y directions. respectively, after the k-th 


iteration: r*( +) 1s the complex conjugate of the data vector; and the step-size parameter 


u, 1S given bv 


fi = y= (Cuenca nO (2.59) 


where uy 1S a constant, and it is equal to the inverse of the sum of the signal and noise 
power, a? and o? , respectively, at the center element of the array. After each iteration 
of the algorithm, each component of the complex weight vectors is normalized by its 
respective magnitude in order to maintain unit magnitude. 

Since We are Working with multiple linear arrays to obtain estimates of the targets’ 
locations, we use average, steady-state, complex weight vectors in the X and Y di- 


rections, that 1s, 
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Ne 


hy (=a), Alan) (2.60) 
n=—N’ 
and 
MM’ 
saya 
by, LD) = 37 2, Bkaem (2.61) 


Note that each of the P x | average, steady-state, complex weight vectors in the X and 


Y directions is given by 


By (4) = [hte (4:1), Pg, (G2) pererrealte, (GP) (2.62) 


avs 


and 


hy (q) =Dhy, (g.1), Bel Deano Cale (2.63) 


Next, if we now use the average, steady-state, complex weight vectors in the X and 
Y directions, then the one-step. steady-state, forward linear prediction errors for har- 


monic number q in the X and Y directions, are given by 


P 
° . ] ° mv my 
F, (4,i,/) = R(q,i,n) — “Pp > hy, (OP)R(QI a ee — — NGO, ee (2.64) 
p=) 


and 
ig 
Fy (qi) = R(g.n,i) — Pp > Fy, Qe) R( qu +p), m=— M,...,0,....40. (2.65) 
p=} 


If we define the forward spatial Z-transform of g(i) as 
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(2.66) 


then 


Zig(i+ p)} = G(z)z ”, (2.67) 


and, upon taking the forward spatial Z-transform of Eqs. (2.64) and (2.65) with respect 
to 1, We obtain 


iE 


l 4 : : 
F, (qin) = R(q.2n)L ~ + Dihx (PEF) nm — NiocccOyn i, ee 
p=1 
and 
i 
F, (qunzy) = R(qmzy)[1— +) Ay, (4 m= — MgO, nee 
Py q; yy ae q, xy P Yeugl — LYE geeegVgereguth « are 
p=! 


From Eqs. (2.68) and (2.69), the transfer functions of the one-stép forward linear 


prediction-error filters (analvsis filters) in the X and Y directions are given by 


F,, (q.2Z,sM) = 
Px 1 7 y 
H pg (4.zx,) = Rigen =]|]- Pp D Prog 4) n=— N’,...,0,...,1) sae 
p= 
and 
F p,(q.m2y) I : 8 ! ' 
Mpg (q.M,2y) = RG) Dp a AqP)z , m=— M....,0,...,47 


Since the average complex weights in the X and Y directions are independent of the in- 


dices n and m, the transfer functions are also independent of n and m, that 1s 


avg 


P 
Hpe (42x) = | -p 2} ong GP )Ex (2.72) 
p=! 


and 


P 


] = 
Hoe (92) =1- Dy, (av)”. 2.73) 
p=l 


Therefore the transfer functions of the inverse filters (synthesis filters) in the X and Y 


directions are given by 


z= : 
A PE(G,2x) —<—<<— (2.74) 
] = 
1) hy (gp )ax” 
p=! 
and 
Hpe,(4.2,) = a. (2.75) 


| —p 
_ Pa ane (q.p )2y 


In order to obtain the spatial-frequency response (angular response) of the filter 
mesembed by Eq. (2./4), let 


+j2rud, 
2, = exp( +)2nfyd,) = exp(———*). (2.76) 
Substituting Eq.(2.76) into Eq.(2.74) yields 
zs 1 
Hor (q,u) = , —— (977) 
] Jempu 
1-- ) Ay,,| (9.p) exp ——— ) 


p= 


Similarly, if we let 


+j2nvady 
zy = exp( tj2af,d, ekp(= =) (2.78) 


in Eq. (2.75), we obtain 


| 


Hpe.(9,¥) = (2.79) 


= —j2npvd, 
pdb Gpal——) 


Equations (2-77) and (2-79) are then evaluated as functions of the direction cosines 
u and v, respectively, for a given harmonic q. The values of u and v that maximize Eqs. 
(2-77) and (2-79), respectively, then become the linear prediction-error estimates u?(q) 
and v?4(q) of u(q) and v,(q) . The estimates of @,(¢) and ~,(qg) at each harmonic q are 
given by [Ref. 3] 





8o(q) = sin ((aPF(qyP + EPA(qyP}?], (2.80) 
and 
‘ to (9) 
j4a)=a0"} — | q#0 (2.81) 
Uy (q) 


where u?§(qg) and r?£(g) are the estimates of the direction cosines obtained by using the 


adaptive linear prediction-error algorithm. 


Il. SIMULATION RESULTS 


Computer simulation results for five test cases are presented in this chapter. The test 


cases Were designed to test both of the algorithms multiple broadband target localization 


capability, full angular coverage capabilities, and their angular resolution as a function 


of the input SNR at a single element in the array, harmonic number, and the number 


of iterations. 


For simplification purposes, in all of the simulation results we use the following 


notation: 


bt 


LMS: When the frequency domain modified complex Least- Mean-Square adaptive 
beamforming algorithm 1s used. 


ALP: When the frequency domain modified complex adaptive linear prediction- 
error filtering beamforming algorithm 1s used. 


The test cases are described below: 


Case | - A single broadband target radiating six harmonics, located at broadside 
relative to the planar array; 


Case 2 - Three broadband targets, each radiating two unique harmonics, located 
at random positions; 


Case 3 - A single broadband target radiating six harmonics, located at endfire rel- 
ative to the planar array; 


Case 4 - Three broadband targets. Two targets are located in the same plane and 
have a common harmonic. Three subcases were examined. Each subcase used a 
different harmonic as the common spectral line. 


. Case 5 - Three broadband targets, located at random positions. Two targets have 


a common harmonic. Three subcases were examined. Each subcase used a different 
harmonic as the common spectral line. 


The simulation results of these five cases are based upon processing the output 


electrical signals from a 11 x 11 planar array of equally spaced hydrophones. The 


acoustic field incident upon the planar array was, in general, the sum of several plane- 


wave fields travelling in different directions. Each plane-wave field consisted of an arbi- 


trary number of harmonics (spectral lines) emanating from one of the broadband sound 
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sources. As a result, the output electrical signal at each element of the planar arrav was 
composed of an arbitrary number Kk of harmonics, all with identical amplitudes of unity. 

The fundamental frequency for all test cases was chosen to be f, = 1000 Hz. Via Eq. 
(2-6), the fundamental period at each element of the array was | millisecond. The sam- 
pling parameter, S, was set equal to 10. The number of time samples taken per element 


of the array 1s given by 


L=SK+5 (3.1) 


where L is the total number of samples taken per element, K is the total number of 
harmonics present in the signal, and S is the sampling parameter. With S= 10 and K=6, 
L=65 time samples per element. 


In the computer simulations presented in this thesis, the interelement spacing 1s 








/ 1 % 
d, = d, = ii , (3.2) 
Where the minimum wavelength 1s 
: a ae 
A 33 
nin = Fa / 
and 
Trae = ez (3.4) 
Substituting Eqs. (3.3) and (3.4) into Eq. (3.2) yields 
a 255 
d, = d, (Kf). (say 


For the ALP algorithm, the order of the forward prediction-error filter was 10. The 
reference element was the first element of each linear array and the number of linear 


arrays that have been processed in both directions was 11. Initial simulation tests indi- 


Zi 


cated that the number of iterations required for the LMS algornthm to approach a 
steady-state error was 100, while for the ALP algorithm 1000 iterations were required. 

Baseline, or “no noise” test case results were the first generated to ensure that both 
algorithms were working properly in a noise-free environment. Figure 3 depicts the 
noise-free, time domain, received signal at element (1,1) for Case 2. For each baseline 
test case, bearing and depression angle estimation errors (measured in degrees) were 
obtained by running the computer simulation once allowing 100 iterations for the LMS 
algorithm and 1000 iterations for the ALP algorithm. 

Following compilation of the baseline results for the first three cases, identical test 
cases Were run using additive, wide-sense stationary, zero mean, white, gaussian noise 
Samples to corrupt the time samples of the received signal. Figure 4 depicts the time 
domain, received signal at element (1.1) for case 2 fora SNR=-9dB. For each test case, 
and for a given input SNR at a single element of the array, average bearing and de- 
pression angle estimation errors were obtained by running the computer simulation 10 
times. During each run. the LMS algorithm was allowed 100 iterations whereas the ALP 


algorithm was allowed 1000 iterations. 
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Figure 3. Real Received Signal at Element (1,1), for Case LMS2 : No Noise 
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Figure 4. Real Received Signal at Element (1,1), for Case LMS2 : SNR=-9dB 
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A. CASE 1 


Case 1 placed a single broadband target at broadside relative to the planar array 
(1.e., 8, = 0°). This was considered the simplest case for both of the algorithms since it 
is at broadside that the far-field beam pattern has its narrowest beamwidth and, as a 
result, both algorithms should provide good angular resolution [Ref 4]. The general 
plane-wave field radiated by the target consisted of six harmonics. 

Figures 5, 6 and 7 present the average estimation errors of the depression angle 
(é,,) versus SNR for Case LMS] at 1000 Hz, 3000 Hz and 6000 Hz, respectively. Results 
for bearing angle estimation errors in the broadside case are irrelevant since the target 
is directly above the array, and the bearing angle has no meaning. From these figures, 
We Can See two trends. Increasing the SNR results in smaller magnitudes of the average 
estimation errors of the depression angle and as the harmonic number(q) increases, the 
magnitude of the estimation error decreases. These two trends were expected. As we 
increase the SNR, the noise components of the received signals are less dominant. An 
increase of the harmonic number (q) represents an increase in the frequency and, as a 
result. the beamwidth of the far-field beam pattern decreases. The decreased beamwidth 
increases the angular resolution of the algorithm and, as result, decreases the estimation 
error. 

Figures 8, 9 and 10 present estimates of the direction cosines U, and V, for Case 
ALPI at SNR=-3 dB and frequencies 1000 Hz, 3000 Hz, and 6000 Hz, respectively. 
Similarly, Figures 11, 12 and 13, and 14, 15, and 16 present estimates of the direction 
cosines U, and V, for Case ALPI with the same order of frequency as above and for 
SNR=-6 dB and SNR=-9 dB, respectively. As we discussed in Chapter 2, the location 
of the peaks of the magnitude of the synthesis filter transfer function in dB give the U, 
and |, estimates. From these figures we can clearly observe that as the SNR decreases, 
the difference between the peak and the background noise in the transfer function plots 
decreases. Also, the harmonic number is important, because increasing the harmonic 
number produces sharper peaks in the transfer function plots resulting in more accurate 
estimation. 

Figures 17, 18 and 19 present the average estimation error of the depression angle 
(€.,) versus SNR for Case ALPI at 1000 Hz, 3000 Hz and 6000 Hz. The trends which 
were present in Case LMS] are apparent here, that is, as SNR increases, the magnitude 
of the estimation error decreases and, as the harmonic number (q) increases, the mag- 


nitude of the estimation error decreases. 
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Finally Tables I and 2 in the Appendix A present the numerical data for Cases 
LMS1 and ALP! respectively for all six harmonics. Comparing the two algorithms us- 
ing these tables, we can conclude that for Case 1, both algorithms give small magnitude 
estimation errors in low SNR with a little better performance for the LMSI particularly 


at 1000 Hz. 
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Figure 18. Average Estimation Error of Depression Angle vs. SNR at 3000 Hz for 


Case ALPI. 
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bee CASE 2 


Case 2 placed three broadband targets at random locations. The general plane-wave 
field radiated by each of the targets contained two unique harmonics . Target | was lo- 
cated at (8, = 49°, w, = 38°) and radiated harmonics | and 6. Target 2 was located at 
sy, = 137°) anderadiated harmonics 2 and 5. Target 3 was located at (0,= 77°, 
W) = 307°) and radiated harmonics 3 and 4. Since there are three incident plane-wave 
fields, each containing two unique harmonics, the output electrical signal from each el- 
ement in the arrav exhibits a total of six harmonics. 

Figures 20 through 25 present the average estimation errors of the depression (é,,) 
and bearing (e,,) angles versus SNR for Case LMS2 for all the harmonics in increasing 
order of the harmonic number (q), respectivelv. The same trend about the SNR as in 
Case | holds here, that 1s increasing the SNR results in smaller magnitude of the average 
estimation error. These general results can be explained using the same arguments pre- 
sented for Case 1. 

Figures 26 through 43 present estimates of the direction cosines U, and V, for all the 
harmonics placed in increasing order of the harmonic number (q) at SNR’s -3 dB, -6 dB 
and -9 dB, respectively. Again, the two trends observed in Case ALPI are apparent here, 
that 1s, decreasing the SNR decreases the difference between the peak and the back- 
ground noise and increasing the harmonic number (q) produces sharper peaks in the 
transfer function plots. 

Figures 44 through 49 present the average estimation errors of the depression (@,,) 
and bearing (e,,) angles versus SNR for Case ALP2 for all the harmonics in increasing 
order of the harmonic number (q), respectively. The same trends as in LMS2 hold, with 
the exception that now the magnitude of the average estimation error for both angles 1s 
smaller than in LMS2. Note the difference in scaling between the error plots for Case 
leevitceand Case ALP2. 

Finally, Tables 3 and 4 in the Appendix A present the numerical data for Cases 
LMS2 and ALP2, respectively, for all six harmonics. Comparing the two algorithms 
using these tables, we can conclude that for Case 2, the ALP algorithm gives smaller 


magnitude estimation errors in most instances than the LMS algorithm especiallv for 
low SNR. 
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at 1000 Hz for Case LM1S2. 
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Figure 23. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 4000 Hz for Case LMS2. 
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Figure 45. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 2000 Hz for Case ALP2. 
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Figure 46. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 3000 Hz for Case ALP2. 
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Figure 47. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 4000 Hz for Case ALP2. 
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Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 6000 Hz for Case ALP2. 
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C. CASES 


Case 3 placed a single broadband target at endfire relative to the planar array (i.e., 
6, = 90°, w, = 90°). This case was considered the most difficult for both algorithms since 
it is at endfire that the far-field beam pattern beamwidth is broadest [Ref. 4]. The general 
plane-wave field radiated by the target consisted of six harmonics. 

Figures 50 and 51 present the average estimation errors of the depression (e,,) and 
bearing (e¢,,) angles versus SNR for Case LMS3 at 1000 Hz and 3000 Hz, respectively. 
The estimation errors for the highest harmonic (6000 Hz) are not shown for Case LMS3 
since the magnitude of these errors were between 45 and 90 degrees. This poor per- 
formance for the highest harmonic at endfire can be explained theoretically [Ref. 6]. 

Figures 52, 53, and 54 present estimates of the direction cosines U, and V, for Case 
ALP3 at SNR=-3 dB and frequencies 1000 Hz, 3000 Hz, and 6000 Hz, respectively. 
Similarly, Figures 55, 56, and 57, and 58, 59, and 60 present estimates of the direction 
cosines U, and V, for Case ALP3 with the same order of frequency as above and 
SNR=-6 dB and SNR=-9 dB. respectively. Here we can observe an ambiguity for the 
}) estimates at 6000 Hz which can be explained by examining the exponential term of 
Eq.(2-79). This term for the highest frequency component creates integer multiplies of 
2z producing two peaks at +1 for the V, estimates at 6000 Hz. 

Figures 61, 62 and 63 present the average estimation errors of the depression (e,,) 
and bearing (e,,) angles versus SNR for Case ALP3 at 1000 Hz, 3000 Hz, and 6000 Hz, 
respectively. We can observe that in Case ALP3 at 6000 Hz we have acceptable esti- 
mates for depression angle (6). In Figure 63 the yw, error plot is not shown since the 
magnitude of this error was between 45 and 90 degrees. 

Finally, Tables 5 and 6 in Appendix A present the numerical data for Cases LMS3 
and ALP3 respectively for all six harmonics. Comparing the two algorithms using these 
tables, we can conclude that for Case 3 the ALP algorithm gives smaller magnitude es- 


timation errors in most instances than the LMS algorithm for both angles. 
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Figure 61. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 1000 Hz for Case ALP3. 
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Figure 62. Average Estimation Errors of Depression aud Bearing Angles vs. SNR 
at 3000 Hz for Case ALP3. 
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Figure 63. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 6000 Hz for Case ALP3. 
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me: CASE 4 


Case 4 corresponds to three broadband targets being present two of which have a 
common harmonic. Target 1 was located at (6, = 45°, Wy, = 0°). Target 2 was located at 
(0, = 45°, y, = 180°). Target 3 was located at ( 6,= 33° , wy, = 47°). A total of six har- 
monics Were present in the output electrical signal. The following three subcases were 
examined: 


e Case 4-A: Three targets, two of which are located in the same plane. Targets | and 
2 have a common spectral line at 1000 Hz. 


e Case 4-B: Three targets, two of which are located in the same plane. Targets | and 
2 have a common spectral line at 3000 Hz. 


e Case 4-C: Three targets two of which are located in the same plane. Targets | and 
2 have a common spectral line at 6000 Hz. 


For these cases, onlv no noise results were compiled for both algorithms in order to 
examine the ability of target discrimination when two targets are in the same plane and 
have a common harmonic. 

Figures 64, 65 and 66 present estimates of the direction cosines U, and V, for each 
common harmonic, that is at 1000 Hz for Case ALP4-A, 3000 Hz for Case ALP4-B and 
at 6000 Hz for Case ALP4-C. Clearly, we can observe that except for case ALP4-A, the 
algorithm can discriminate the location of two targets. 

Tables 7 through 12 in the Appendix A present numerical data for all the subcases 
of Case 4 for both algorithms. From these tables we can observe that for all the har- 
monics, except the common one, both algorithms correctly identified the location of the 
targets with zero degrees estimation errror. However, for the common harmonic, the 
LMS algorithm located a false fourth target, which was located exactly between the ac- 
tual locations of targets 1 and 2. On the other hand the ALP algorithm, except for the 
subcase ALP4-A, for the common harmonic, identifies two different estimates of the di- 
rection cosine U, from the transfer function plots, which results in two different target 


locations. 
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Figure 66. 


me CASE 5 


Case 5 placed three broadband targets at random locations. Target 1 was located 
at (8, = 49°, wy, = 38°). Target 2 was located at (@, = 5°, yw, = 137°). Target 3 was located 
at (0, = 77°, w, = 307°). A total of six harmonics were present in the output electrical 
signal. The following three subcases were examined: 


e Case 5-A: Three targets located at random positions. Targets 1 and 2 have a 
common spectral line at 1000 Hz. 


e Case 5-B: Three targets located at random positions. Targets 1 and 2 share a 
common spectral line at 3000 Hz. 


e Case 5-C: Three targets located at random positions. Targets 1 and 2 share a 
common spectral line at 6000 Hz. 


For these cases, onlv no noise results were compiled for both algorithms in order to 
examine the ability of target discrimination when two targets are in random positions 
and have a common harmonic. 

Figures 67, 68 and 69 present estimates of the direction cosines U, and V, for each 
common harmonic, that is. at 1000 Hz for Case ALP5-A, at 3000 Hz for Case ALP5-B 
and at 6000 Hz for case ALP5-C. Clearly, we can observe that exept for case ALPS5-A, 
the algorithm can discriminate the location of two targets. 

Tables 13 through 18 in the Appendix A present numerical data for all the subcases 
of Case 5 for both algorithms. From these tables we can observe that for all the har- 
monics, except the common one, both algorithms correctly identified the location of the 
targets with zero degrees estimation errror. However, for the common harmonic, the 
LMS algorithm located a false fourth target. On the other hand the ALP algorithm, 
except for subcase ALPS-A, for the common harmonic, identifies two different estimates 
of the direction cosines U, and V, from the transfer function plots, which results in two 


different target locations. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The goal of this thesis was to evaluate the multiple broadband target localization 
capabilitv and the full angular coverage capability of two frequency domain adaptive 
beamforming algorithms. The main advantage of both algorithms is that no a priori 
knowledge is assumed about the statistics of the received signal, the total number of 
targets present, where they are located, and the frequency content of the targets. 

The two frequency domain adaptive beamforming algorithms evaluated were the 
modified complex Least-Mean-Square (LMS) algorithm developed by Ziomek and Chan 
(Ref. 6], and the modified complex Adaptive Linear Prediction-Error (ALP) algorithm 
which was presented for first trme. Chapter 3 presented the results of five test cases 
which were designed to test both algorithms for several different capabilities. It should 
also be noted here that the number of time samples taken per element of the arrav never 
exceeded 65 time samples, an amount which represents a small number of data points. 

Summarizing the simulation results discussed in Chapter 3, we can list the following 


conclusions: 


e In most cases tested, both algorithms can localize multiple broadband targets and 
they have full angular coverage capability. 


e Increasing either the sampling rate or the array size vields better estimates in lower 
SNR for both algorithms 


¢ ALP generally gives better estimates in lower SNR than LMS 


¢ ALP can, in most of the “no noise” cases, discriminate between two targets that 
have a common harmonic while LMS cannot 


¢ ALP is more computationally intensive than LMS and requires much more iter- 
ations to approach a steady-state error. 


¢ ALP is more flexible than LMS because it can process only a part of the total 
number of array elements to localize the target 


In the course of this investigation, several possible areas for future research pre- 


sented themselves: 


OW 


further study of the effects of varving the step-size parameter, uw , in an attempt to 
decrease the magnitude of the estimation error. 


investigation of parallel processing implementation for the ALP algorithm in order 
to reduce the computation time. 


investigation of larger array size in lower SNR mainly using the ALP algorithm 


investigation of a noise reduction system prior to processing the output electrical 
signals. 
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APPENDINX A. 


This appendix presents the numerical data for all of the test cases. Tables | through 
6 present numerical data for the three first “noisy” cases, Cases 1 through 3, for both 
algorithms, while tables 7 through 12 present numerical data from the two “no noise” 
cases, Case 4 and 5, for both algorithms. The following notation is used in all the tables: 


e LMS: When the frequency domain modified complex Least- Mean-Square adaptive 
beamforming algorithm is used. 


ALP: When the frequency domain modified complex adaptive linear prediction- 
error filter beamforming algorithm 1s used 


e §@,: Depression Angle 
e w,: Bearing Angle 
e é,,: Average Estimation Error for Depression Angle 


e ¢,,: Average Estimation Error for Bearing Angle 
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Table 1. NUMERICAL DATA CORRESPONDING TO CASE LMSI1 


Case 1: One target located at broadside relative to planar array - six harmonics 
present in output electrical signal - LMS algorithm 1s used - 100 iterations 


Location (deg) Average Estimation Error (deg) 
ms a= zm 
nioni 
aa a 


ts a 
a OO 
[sot 0.561 -0.607 12383 
[0 | 270 | 0360 | 0.760 | 0.689 
[si to 029 0.596 
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Table 2. NUMERICAL DATA CORRESPONDING TO CASE ALP! 


Case 1: One target located at broadside relative to planar array - six harmonics 
present in output electrical signal - ALP algorithm is used - 1000 iterations 


SXR R 
—— 
a a 


CO os 
2 ato 307 2.392.308 
0 ts st 1.993 
A 
OS 
et oar 0.665 0981 
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Table 3. NUMERICAL DATA CORRESPONDING TO CASE LMS2 


Case 2: Three targets located at random positions - six harmonics present in output 
electrical signal - LMS algorithm 1s used - 100 iterations 


Location (deg) Average Estimation Error (deg) 
Harmonic 
(q) 


a eT 
Pn | tn fe | |e | a 


a | 26s [oa [oan [0832 [os [0335 
[2s | a7 _| 0077 [27a [oer [-s017 | 0097 [3.051 
C3 tt] 307_| 3.037 | -0.108 | 0351 [2.025 [7991 [9333 
aot aati cert Lesose [oni [ase [sense 
Te sf .0s assy [10.857] 4382 [2355 [8.780 
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Table 4. NUMERICAL DATA CORRESPONDING TO CASE ALP2 


Case 2: Three targets located at random positions - six harmonics present in output 
electrical signal - ALP algorithm 1s used - 1000 iterations 


| Location (deg) | | Location (deg) | Average Estimation Error (deg) 
Harmonic SS 
) 
oe 


sate [00s [0.896 [0380 [283 2a 
a EE 
3307 fires [00 [3.207 [0361 [0.697 1st | 
a at_ 307.968 [0.526 [0.964 [0.218 [1.763 [0.035 
PS sis Poo anor [one 25t [o50 f-5 91 

[asf 38 [0.173 [0.146 [-0.197 | -0.019 | 2.068 [5.109 
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Table 5. NUMERICAL DATA CORRESPONDING TO CASE LMS3 


Case 2: One target located at endfire relative to planar array - six harmonics present 
in output electrical signal - LMS algorithm 1s used - 100 iterations 


— 

) Cb 
71.186 
[2507901 33.033 | -0.200 | 40.295] -0.646 | 27.236 | -0330_ 
[390 [9019.90 J -0.019 | 16.401 | -0.229 | 16.706 | -0.325. 
[90 902.257 [0.006 | 7.365 | 0.287 | 12.025 | 0.139 
Ts [990 2.136 | 0083 [1.146 | 0.098 | 15.699 [18.528 
[6 [90 [9090877 85.120) 69.826 | 64817] 46.737 [90.393 
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Table 6. NUMERICAL DATA CORRESPONDING TO CASE ALP3 


Case 2: One target located at endfire relative to planar array - six harmonics present 
in output electrical signal - ALP algorithm is used - 1000 iterations 


Location (deg) Average Estimation Error (deg) 


I, 


905948 | 0217 | 806 |0136 | oor? | 098 
[290 [907.287 [0.362 [11.323 0.329 

[3 [90 [903342 | -0.107 | 3.732 |-0370 | 1.935 [0.767 
Cf 90 | 902.233 [os | 21as [0286 [7.725 [0.025 
[3 [90] 90] 0.986 [-0.160 | 3.512 | 0076 [14.090 | 6512 
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Table 7. NUMERICAL DATA CORRESPONDING TO CASE LMS4-A 


Case da: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 1000 Hz - six harmonics present in output elec- 
trical signal - LMS algorithm is used - 100 iterations - no noise 


Location (deg) Average eee oe Error 
Harmonic (q) (deg) 


ee ee, 
[as 0 | 45.000 90.000 

13 
i 
6s 1 a | 090 | cm 
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Table 8. NUMERICAL DATA CORRESPONDING TO CASE ALP4-A 


Case 4a: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 1000 Hz - six harmonics present in output elec- 
trical signal - ALP algorithm 1s used - 1000 iterations - no noise 


: Location (deg) Average Bo oauOn EmnoG 
Harmonic (q) (deg) 


7 

[0 [sare 
[180 [2.921 |_180.000 
as 0008 00 

Cat | 0037 [0.011 
Tat 0.037 001 
ores at 00 ot 













U2 
a eh Ud a = a 
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Table 9. NUMERICAL DATA CORRESPONDING TO CASE LMS4-B 


Case 4b: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 3000 Hz - six harmonics present in output elec- 
trical signal - LMS algorithm 1s used - 100 iterations - no noise 


Harmonic (q) (deg) 


ae 
i 
[745.00 [90.00 
ists son 
[at | 0.000 | 0.000 
Ee 
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Table 10. NUMERICAL DATA CORRESPONDING TO CASE ALP4-B 


Case 4b: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 3000 Hz - six harmonics present in output elec- 
trical signal - ALP algorithm is used - 1000 iterations - no noise 


| 
[tame | 
a 
0 | 0009 [0.000 
so se to 
iso 156 [0.00 
Se os 
180 [0.009 [0.000 
at nt or 


















Cd td (2 
td a oe) a i a Lad 
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Table 11. NUMERICAL DATA CORRESPONDING TO CASE LMS4-C 


Case 4c: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 6000 Hz - six harmonics present in output elec- 
trical signal - LMS algorithm is used - 100 iterations - no noise 


Location (deg) Average Estimation Error 
Harmonic (q) (deg) 


ee 


=a [000] ooo 
[0 | 48.000 [90.000 
| a5] 180 | aso —t sam 


SS 
Gr a 
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Table 12. NUMERICAL DATA CORRESPONDING TO CASE ALP4-C 


Case 4c: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 6000 Hz - six harmonics present in output elec- 
trical signal - ALP algorithm 1s used - 1000 iterations - no noise 


. Location (deg) Average San Error 
Harmonic (q) (deg) 


a 
as 
[as [180 0.009 [0.000 
Mess «| 0 ef 2070-4 0.000 
[SSCS~sCSS 802.070 [0.000 
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Table 13. NUMERICAL DATA CORRESPONDING TO CASE LMS5-A 


Case 5a: Three targets located at random positions - targets one and two have a 
common spectral line at 1000 Hz - six harmonics present in output electrical signal - 
LMS algorithm is used - 100 iterations - no noise 


Average Estimation Error 
Location (de 
as 


a 
a 
[7307 0.000 [0.00 
aso 0000 
[3137 | 0.00] 0.000 
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Table 14. NUMERICAL DATA CORRESPONDING TO CASE ALP5-A 


Case 5a: Three targets located at random positions - targets one and two have a 
common spectral line at 1000 Hz - six harmonics present in output electrical signal - 
ALP algorithm is used - 1000 iterations - no noise 












Average Estimation Error 










Location (deg) 
Harmonic (q) 


oo 
ss 926 
a a 
[3st P0006 | 0.328 





tt 
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Table 15. NUMERICAL DATA CORRESPONDING TO CASE LMS5-B 


Case 5b: Three targets located at random positions - targets one and two have a 
common spectral line at 3000 Hz - six harmonics present in output electrical signal - 
LMS algorithm is used - 100 iterations - no noise 


Average Estimation Error 
Location (de & 
Harmonic (q) (deg) 
a a nn, lee, | omc | 


[| 3070.00 [0.000 
a 
ss i | at 000 | ten 
Taco 000 0.000 
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Table 16. NUMERICAL DATA CORRESPONDING TO CASE ALP5-B 


Case 5b: Three targets located at random positions - 
targets one and two have a common spectral line at 3000 Hz - six harmonics present 
in Output electrical signal - ALP algorithm 1s used - 1000 iterations - no noise 


Location (deg) Estimation Error Se 
peers aa 
es 
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Table 17. NUMERICAL DATA CORRESPONDING TO CASE LMS5-C 


Case 5c: Three targets located at random positions - targets one and two have a 
coniumon spectral line at 6000 Hz - six harmonics present in output electrical signal - 
LMS algorithm is used - 100 iterations - no noise 


; Average Estimation Error 
Location (deg 
Harmonic (q) ( s) < _— 




















[sf st [0.000 [0.000 
st 90 
15.000 | 252.394 


Table 18. NUMERICAL DATA CORRESPONDING TO CASE ALP5-C 


Case 5c: Three targets located at random positions - targets one and two have a 
common eeu line at 6000 Hz - six harmonics present in output electrical signal - 
ALP algorithm 1s used - 1000 iterations - no noise 


. - g S . * 
Location (deg) Average E ree Error 
Harmonic (q) (deg) 


a 
[77 07 | 0.098 | 0.012 
a ee ee 
[ss SC~S ti a «da es 
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